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Abstract 



Sleptons, neutralinos and charginos were searched for in the context of scenar- 
ios where the lightest supersymmetric particle is the gravitino. It was assumed 
that the stau is the next-to-lightest supersymmetric particle. Data collected 
with the DELPHI detector at a centre-of-mass energy near 189 GeV were anal- 
ysed combining the methods developed in previous searches at lower energies. 
No evidence for the production of these supersymmetric particles was found. 
Hence, limits were derived at 95% confidence level. 
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1 Introduction 

Supersymmetry (SUSY) may be broken at a scale below the grand-unification scale 
Mgut, with the ordinary gauge interactions acting as the messengers of supersymmetry 
breaking |]T],0. In the corresponding models (GMSB models), the gravitino, G, turns out 
to be the lightest supersymmetric particle (LSP) and is expected to be almost massless. 
The next-to-hghtest supersymmetric particle (NLSP) is therefore unstable and decays, 
under the assumption of i?-parity conservation, into its ordinary matter partner and an 
invisible gravitino. 

The number of generations of supersymmetry breaking messengers and the value of 
tan/? usually determine which supersymmetric particle is the NLSP 0,^J^J^. In the 
majority of the GMSB space, the NLSP is a slepton, /. Moreover, depending on magnitude 
of the mixing between tr and f^, there exist two possible scenarios. If the mixing is 
large [|, fi is the NLSP, but if the mixing is negligible, fi is mainly right-handed and 
almost mass degenerate with t_he other sleptons. In this case, the and jln three body 
decay (/-^ fir/ with r G), is very suppressed, and cr and jln decay directly into 

IG. This scenario is called / co-NLSP. Searches for supersymmetric particles within both 
these scenarios are described in this article. 

Due to the coupling of the NLSP to G, its mean decay length can range from fim 
to meters depending on the mass of the gravitino {ttiq), or equivalently, on the scale of 
SUSY breaking, y/P: 




1.76 X 10-3,/ ( ^ ) - 1 ( ) ' (^)' cm, (1) 

For example, for ttiq < 250 eV (a/F < 1000 TeV), the decay of the NLSP can take 
place within the detector. This range of \/F is in fact consistent with astrophysical and 
cosmological considerations 

The results of three searches are presented in this work. The first one looks for the 
production of Xi pahs with either Xi decaying to fir and fi then decaying promptly 
into tG, which is an update of the search presented in ref. |]10|, or Xi decaying to // 



with BR(x5 ^0 = 1/3 and I promptly decaying into IG: e+e^ XiXi ~^ ~^ 
IGll'GV . These two modes represent the two extremes in the range of possible decays of 
the neutralino. In particular, a higgsino-like x? would decay only to fir for all practical 
purposes since the higgsino component of the x? couples to / through Yukawa couplings. 
On the other hand, the decays of a gaugino-like Xi are regulated only by phase space 
considerations. Therefore, in the case of the fi-NLSP scenario, neutralino pair production 
would mainly lead to a final state with four tau leptons and two gravitinos, while in the 
case of a co-NLSP scenario, the final signature would contain two pairs of leptons with 
possibly different flavour and two gravitinos. 

The second search concerns / pair production followed by the decays / IG within the 
detector volume. The signature of such an event will be at least one track of a charged 
particle with a kink or a decay vertex when the / decays inside the tracking devices. If 
the decay length is too short (small m^) to allow for the reconstruction of the I track, 
only the corresponding lepton or its decay products will be seen in the detector, and the 
search will then be based on track impact parameter. However, if the decay takes place 
outside the tracking devices (large m.^), the signature will be that of a heavy charged 

^In GMSB models large mixing occurs generally in regions of tan/3 > 10 or > 500 GeV. 
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particle already studied by DELPHI [Ol. For very light gravitinos the limits from the 



search for sleptons in gravity mediated (MSUGRA) models can be applied [12 . 

The third search looks for the pair-production of lightest charginos, xf. Charginos, 



if produced in this context, would promptly decay through the channel Xi ^i^^t |T3 
The fi would then decay into rG, with non negligible mean lifetime. This search is 
divided into three sub-channels according to the mean lifetime of the stau as explained 
in the previous paragraph: two acoplanar leptons with missing energy, at least one track 
with large impact parameter or a kink, or at least one track corresponding to a very 
massive stable particle. 

The data samples and event selections are respectively described in sections |^ and ^, 
while the results and a model dependent interpretation are presented in section 



2 Event sample and experimental procedure 

All searches are based on data collected with the DELPHI detector during 1998 at a 
centre-of-mass energy of 189 GeV. The total integrated luminosity was 153.6 pb-\ A 
detailed description of the DELPHI detector can be found in and its performance in 



15| . In all cases, the ri-NLSP scenario searches are updates to similar searches carried 
out at lower centre-of-mass energies. All co-NLSP scenario searches are carried out at 
= 189 GeV. 

To evaluate the signal efficiencies and Standard Model (SM) background contami- 



nations, events were generated using different programs, all relying on JETSET 7.4 |I6 



tuned to LEP 1 data JT^ for quark fragmentation. The program SUSYGEN ||T8| was used 



to generate the neutralino pair events and their subsequent decay products. In order 
to compute detection efficiencies, a total of 42000 events were generated with masses 
67 GeV/c2< rur^ + 2 GeV/c^ < m^o < ^/s/2. A f pair sample of 36000 events (subdi- 
vided in 36 samples) was produced with PYTHIA 5.7 |TB| with staus having mean decay 



lengths from 0.25 to 200 cm and masses from to 90 GeV/c^. Another sample of f 
pair was produced with SUSYGEN for the small impact parameter search with irif^ from 
7 to 80 GeV/c^. Similar samples of smuons and selectrons were produced to study the 
sleptons co-NLSP scenario. 

SUSYGEN was also used to generate the xf pair production samples and their decays. 
In order to compute detection efficiencies, a total of 45 samples of 500 events each were 
generated with G masses of 1, 100 and 1000 eV/c^, mf-^ +0.3 GeV/c^ < m-+ < \/s/2 and 
> 65 GeV/c^. Samples with smaller Am = m-+ — nif^ were not generated because 

in this region the xf does not decay mainly to fi and Vr but into W and G. The different 
background samples and event selections are described in references [|l9| , ^pA| for rriQ =1 
and 1000 eV/c^ respectively. For the case of uiq = 100 eV/c^, the analysis is the same as 
the one used in the search for sleptons of this paper and consequently, the same sample 
of backgrounds is used. 

The background process e"'"e~ QQin'j) was generated with PYTHIA 5 . 7, while DYMU3 



2T| and KORALZ |2^ were used for /i+/i~(7) and r"'"r~(7), respectively. The generator of 
reference I^Bj was used for e+e~ e^e~ events. 

Processes leading to four-fermion final states, {Z / •y) {Z / •y) , W+W~, Wei^e and Ze+e~, 
were also generated using PYTHIA. The calculation of the four-fermion background was 



verified using the program EXCALIBUR [Q, which consistently takes into account all 
amplitudes leading to a given four-fermion final state. 
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Two-photon interactions leading to hadronic final states were generated using 

QPM and QCD components. The generators of 



TWOGAM 25 



separating the VDM 
Berends, Daverveldt and Kleiss ||26| were used for the leptonic final states. 

The cosmic radiation background was studied using the data collected before the 
beginning of the 1998 LEP run. 

The generated signal and background events were passed through the detailed simu- 
lation |]T3| of the DELPHI detector and then processed with the same reconstruction and 
analysis programs used for real data. 



3 Data selection 

3.1 Neutralino pair production 



The selection used in the search for the process e~^e XiXi ~^ tittit tGttGt 



has been described in |T0|. A very similar selection was used in the general search for 
e+e" ^ XiXi IIVV IGll'GV within the co-NLSP scenario, where BR(x? ^ //) = 
1/3 for each leptonic flavour. The main two differences between these two cases comes 
from the fact that the mean number of neutrinos carrying away undetected energy and 
momentum and the number of charged tracks per event is considerably bigger for the 
fi-NLSP scenario. These differences can be appreciated in fig. |I|, where the simulated 
missing energy normalized to the centre-of-mass energy and the number of charged tracks 
are represented for events with neutralinos weighing 82 GeV/c^, and staus of 80 GeV/c^, 
compared to events with same-mass neutralinos and degenerate sleptons of 80 GeV/c^. 

The pre-selection of events is common to both scenarios, and has been described in |1lO| , 
together with the selection of the search for x\ ~^ ^I'^'i which has not been changed for 
the present analysis. Only the details of the search II are presented in the following. 
Two sets of cuts were applied in order to reduce the 77 and ff(7) backgrounds and a 
third set of cuts to select events according to their topology: 

1- Cuts against 77 backgrounds: the transverse energy, Et, should be bigger than 4 
GeV. The energy in a cone of 30° around the beam axis was further restricted to 
be less than 60% of the total visible energy to avoid possible bias from the Monte 
Carlo samples. The missing mass should be smaller than 0.88a/s. The momentum 
of the charged particle with largest momentum should be bigger than 8 GeV/c. 
The transverse missing momentum, px, should be bigger than 6 GeV/c. These cuts 
reduced the 77 background by a factor of the order of 40. 

2- Cuts against ff (7) and 4-fermion backgrounds: the number of good tracks should be 
smaller than 7. The maximum thrust was further reduced from 0.99 to 0.95. Dividing 
each event into two jets with the Durham algorithm, its acoplanarity should be bigger 
than 8°. The missing mass of the events should be bigger than 0.2-^5. After these 
cuts, the ff (7) and 4-fermion backgrounds were reduced by a factor of the order of 
30. 

3- Cuts based on topology: signal events tend naturally to cluster into a 4-jet topology. 
When events are forced into a 4-jet configuration, all jets should be at least 18° 
away from the beam direction. When reduced by the jet algorithm into a 2-jet 
configuration, the charged particles belonging to each of these jets should be in a 
cone broader than 25°. Finally, the axes of each of the four jets should be separated 
from the others at least by 9°. 
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Figures to ^ show some of the distributions relevant for these selection criteria. The 
discrepancy between data and simulation on the last two bins of figure ^a is in a region 
of soft 77 events that is not relevant to the final results. Tables |1| and ^ show the effect 
of these cuts on the data, expected background and the signal for m^o = 87 GeV/c^ and 

nir^ = 75 GeV/c^ for the cases Xi ~^ Xi ^ 't^ respectively. 

One event was observed to pass the search for neutralino pair production in the fi- 
NLSP scenario, with 1.16 ± 0.19 SM background events expected. Two events pass 
the search for neutralino pair production in the co-NLSP scenario, with 1.2 ± 0.30 SM 
background events expected. After these cuts, efficiencies between 27.0 and 40.7% were 
obtained for the signal events. 



3.2 Slepton pair production 

This section describes the update of the search for the process e^e~ fifi tGtG 
already described in |ll0| , p7| . An additional 153.6 pb~^ integrated luminosity collected 



at the centre-of-mass energy of 189 GeV has been analysed using the same procedure as 
for the data collected at 183 GeV and using the same values for the data selection cuts. 
The same selection cuts have been apphed to the search for cr- and fiR- pair production 
in the framework of / co-NLSP scenario. Therefore, only results and efficiencies will be 
reported in this section, since the details of the selection criteria can be found in pU], 



3.2.1 Search for secondary vertices 



This analysis exploits a peculiarity of the / —>■ IG topology in the case of intermediate 
gravitino masses (i.e. 0.5 eV/c^ < < 200 eV/c^ as dictated by eq. P, namely, one 
or two tracks coming from the interaction point and at least one of them with either a 
secondary vertex or a kink. 

Rather loose preselection cuts were imposed on the events in order to suppress the low 
energy background (beam-gas, beam- wall, etc), 77, e^e~ and hadronic events. The only 
cut that was changed with respect to the analysis at 183 GeV is the total electromagnetic 
energy required in the event. It was increased to ^/i/2 in order to improve the efficiency 
for selectrons. This did not increase noticeably the background contamination by Bhabha 
events. These preselection cuts left about 0.6% of the whole data sample. The events 
that survived the preselection cuts underwent the search for secondary vertices or kinks. 

Fake decay vertices could be present amongst the reconstructed secondary vertices, 
being produced by particles interacting in the detector material or by radiated photons 
if the particle trajectory was reconstructed into two separated tracks. To eliminate these 
classes of events, additional requirements were imposed: 

- to reject hadronic interactions, any reconstructed hadronic interaction (secondary 
vertices reconstructed in region where there is material) must be outside a cone of 
half angle 5° around the slepton direction; 

- to reject segmented tracks, the angle between the tracks used to define a vertex had 
to be larger than 6°; 

- to reject photon radiation in the case of r clusters with only one track, there had 
to be no neutral particle in a 3° cone around the direction defined by the difference 
between the fi momentum and the momentum of the r daughter calculated at the 
crossing point. 
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If no pair of tracks was found to survive these conditions, the event was rejected. Fi- 
gure ^ shows the distribution of these three quantities. The distributions compare real 
data, expected Standard Model background simulation and simulated signal for m^^ = 
60 GeV/c^ decaying with a mean decay length of 50 cm. The excess of data in the first 
bins of fig. |^b is due to an underestimation in the simulation of mismatchings between 
the tracking devices. 

One event in real data was found to satisfy all the conditions described above, while 
l.lSl^Qgg were expected from SM backgrounds. The event was compatible with a 77 — 
r+r^ with a hadronic interaction in the ID detector. 

The vertex reconstruction procedure was sensitive to radial decay lengths, R, between 
20 cm and 90 cm. Within this region a vertex was reconstructed with an efficiency of 
~52%. The VD (Vertex Detector) and the ID (Inner Detector) were needed to reconstruct 
the fi track and the TPC (Time Projection Chamber) to reconstruct the decay products. 
The shape of the efficiency distribution was essentially flat as a function of R, going 
down when the fi decayed near the outer surface of the TPC, due to inefficiencies in the 
reconstruction of the tracks coming from the desintegration products of the r. Also, the 
sensitive region and the efficiency of the vertex reconstruction at 189 GeV was slightly 
lower than at 183 GeV due to the loss of tracks not reconstructed when their VD hits had 
no information in the z direction. Such tracks were not reconstructed for the 189 GeV 
run. However, some of the efficiency lost in the vertex search was recovered later by the 
search based on large impact parameter. 

The search for vertices had an efficiency of the order of 46% for f\ masses between 40 
and 85 GeV/c^ with a mean decay length of 50 cm. The efficiencies decreased near the 
kinematical limit due to a small boost that allowed for big angles to appear between the 
T\ and the desintegration products of the r. For fi masses below 40 GeV/c^, the efficiency 
decreased gradually due to the cut that rejects segmented tracks. This happened because 
the resulting big boost causes the angle between f\ and r decay products to be very small. 

As already said, the same selection criteria was applied to smuons and selectrons. The 
efficiency for selectrons decreased due to the preselection cut on total electromagnetic 
energy (lower than 0.5-ys) and was around 31% for vn^^ between 40 to 85 GeV/c^, while 
the smuon efficiency increased to 55% for the same mass range. 

3.2.2 Large impact parameter search 



To investigate a region of lower gravitino masses the previous search was extended to 
the case of sleptons with mean decay length between 0.25 cm and approximately 10 cm. 
In this case the / track is not reconstructed and only the / (or the decay products in 
the case of f) is detected. The impact parameter search was only applied to those events 
accepted by the same general requirements as in the search for secondary vertices, and not 



selected by the vertex analysis. The same selection criteria described in references ||10| , |27 
were apphed. 

The efficiencies were derived for the different fi masses and decay lengths by applying 
the same selection to the simulated signal events. The maximum efficiency was 32% 
corresponding to a mean decay length of 2.5 cm. The efficiency decreased very fast for 
lower decay lengths due to the cut on minimum impact parameter. The efficiency at 189 
GeV was slightly larger than at 183 GeV since some events not passing the secondary 
vertex selection were recovered in this search, as explained before. For longer decay 
lengths, the appearance of reconstructed f\ tracks in combination with the cut on the 
maximum amount of charged particle tracks caused the efficiency to decrease smoothly. 
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This decrease is compensated by a rising efficiency in the search for secondary vertices. 
For masses above 30 GeV/c^ no dependence on the fi mass was found far from the 
kinematic hmit. However for lower masses the efficiency decreased and it was almost zero 
for a 5 GeV/c^ fi. 

The same selection was applied to selectrons and smuons. For smuons the efficiency 
increased to 59% for a mean decay length of 2.5 cm and masses over 30 GeV/c^ since the 
smuon has always one prong decay. For selectrons the efficiencies were almost the same 
as those for staus. 

Trigger efficiencies were studied simulating the DELPHI trigger response to the events 
selected by the vertex search and by the large impact parameter analysis, and were found 
to be around 99%. 

No events in the real data sample were selected with the above criteria, while 0.321q'Jq 
were expected from SM backgrounds. The number of expected background events at 
\fs = 189 GeV is shown in Table ^ for the combination of the vertex and large impact 
parameter searches. 

3.2.3 Small impact parameter search 



The large impact parameter search can be extended further down to mean decay 
lengths of around 0.1 cm. The same selection criteria described in reference was 
applied. However, some extra selection was added in order to reduce background from 
detector noise or failure, cosmic radiation and rr events. 

Events with anomalous noise in the TPC were rejected requiring less than 20 charged 
particles (before track selection) and relative error of the measured momentum of the 
leading tracks (charged particles with largest momentum in each hemisphere) less than 
50%. The cosmic muon rejection was improved by requiring that the leading tracks with 
impact parameters larger than 1 cm must be reconstructed in the TPC. To reduce the rr 
background where one of the taus decays into a three prong topology (when the single 
track is not reconstructed ), and also gamma conversions, any leading track must have 
at least other charged particle at an angular distance larger than 5°. 

The efficiency of the search turned out not to depend on the fi mass for masses over 
40 GeV/c^, but rather on the fi decay length in the laboratory system. The maximum 
efficiency was ~ 38% for a mean decay length of ~ 2 cm, the efficiency dropped at small 
decay lengths (~ 15% at 1 mm). 

The same selection criteria were used to search for smuons. To search for selectrons, 
in order to increase efficiency, the cut {Ei + E2) < 0.7Ebeam (where Ei, E2 are the 
electromagnetic energy deposits associated to the leading tracks) was not applied. The 
Bhabha events that survived the selection, when the previous rejection cut was not ap- 
plied, were those where at least one of the electrons underwent a secondary interaction, 
thus acquiring a large impact parameter. However, it was found that in these cases the 
measured momentum of the electron was smaller than the electromagnetic energy depo- 
sition around the electron track. Therefore, the cut {Ei/pi + E2/P2) < 2.2 was used for 
the selectron search. The maximum efficiency reach for the smuon search was 43% and 
for the selectron search 35% at 2 cm mean decay length. 

The number of events selected in data was 4, and 4:.54:~^qIj events were expected from 
Standard Model background (see table Two of the candidates had tracks with ffiting 
problems and the other two events were compatible with Standard Model rr events. 
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3.3 Chargino pair production 

The analyses used to search for the hghtest charginos varies according to the stau 
hfetime. Figures | and ^ illustrate the distributions of some of the main variables used 
in the analyses described respectively in section (the stau decays with big impact 
parameter or producing a kink), and |jl9| (the stau decays at the main vertex). The 
plots show real data, expected Standard Model background, and a simulated signal of 
m-+ = 85 GeV/c^ and m-+ = 69 GeV/c^. For the three analyses described in section ||, 
Table ^ shows the range of efficiencies, the main components and the total amount of the 
expected background events, and the number of observed data events for each sample. 



4 Results and interpretation 

Since no evidence for a signal was found in the data, limits on the cross-section of spar- 
ticle pair production were derived. In what follows, the model described in reference |Q 
will be used in order to derive limits. This is a model which assumes radiatively broken 
electroweak symmetry and null trilinear couplings at the messenger scale. The SUSY 
soft parametes, gauge and Yukawa couplings are evolved between the electroweak scale 
(chosen to be rrit) and the messenger scale following the prescription of |2^. The masses 



of gauginos and sfermions at the messenger scale are calculated taking into account cor- 
rections arising from threshold effects. The corresponding parameter space was scanned 
as follows: 1 < n < 4, 5 TeV < A < 90 TeV, 1.1 < M/A < 10^ 1.1 < tan/5 < 50, and 
sign{^) = ±1, where n is the number of messenger generations in the model, A is the 
ratio between the vacuum expectation values of the auxiliary component and the scalar 
component of the superfield and M is the messenger mass scale. The parameters tan j3 
and n are defined as for MSUGRA. 

4.1 Neutralino pair production 

Limits on the cross-section for neutralino pair production were derived in the two 
scenarios for each (m^o,m|-) combination. For the fi-NLSP case, the combination took 
also into account the results from the LEP runs of 1996 and 1997 []10| . 



Figure 0-a shows the 95% C.L. upper limit for the Xi P^ir production cross-section at 
^/s = 189 GeV as a function of m^o and m^^ after combining the results of the searches 



from ^/s = 161 up to 189 GeV with the likelihood ratio method [29|, and scanning 
through the whole parameter space. Figure |^b shows the 95% C.L. upper limit for the 
Xi pair production cross-section at ^/s = 189 GeV as a function of m^o and ttlj^. For 
different number of messenger generations, the ratios between production cross sections 
at different energies are bound to vary within certain limits. Figure ^ presents as an 
example the case of n = 3. For the other scenarios considered in this study (n = 1, 2 
and 4), the maximum difference with respect to figure |^ is not bigger than 10%. This 
variation 



4.2 Slepton pair production 

Figure H shows the 95% C.L. upper limit on the slepton pair production cross-section at 
^/s = 189 GeV after combining the results of the searches at ^/s = 130-189 GeV with the 
likelihood ratio method p9|. The results are presented in the {irtQ^mi) plane combining 
the two impact parameter searches, the secondary vertex analysis and the stable heavy 
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lepton search [|r^ . In particular, figure g-a shows that the minimum upper hmits achieved 
for a given fi were around 0.05-0.10 pb depending on itlq. For ttiq > 9 eV/c^ and a 
80 GeV/c^ fi, a 0.10 pb hmit was obtained. Figures ||-b and -c show the corresponding 
upper hmit for fiR- and pair production cross-sections. Assuming mass degeneracy of 
the three supersymmetric particles, fi, cr and fiR, figure ||-d shows the 95% C.L. upper 
limit for the Ir pair production cross-section. 

4.3 Chargino pair production 

Limits on the production cross-section for chargino pairs were derived for each 
(m(j,mf^,m-+) combination. Figure ^ shows the 95% C.L. upper limit on the chargino 
pair production cross-section at y/s = 189 GeV as a function of m-+ and rrif^ after com- 
bining the results of the searches at ^/s = 183 and 189 GeV with the maximum likelihood 
ratio method |2| for Am > 0.3 GeV/c^ and = 1, 100 and 1000 eV/c^. These lim- 



its, which directly reflect the efficiencies of the applied selections, can be understood as 
follows: 

ttlq = 1 eV/c^: The efficiency of this analysis depends mainly on the mass of the 

chargino. To smaller chargino masses correspond bigger event 
missing energies, and bigger efficiencies. 

ttiq = 100 eV/c^: The map of efficiencies is the result of the convolution of two 

factors. First, larger stau masses imply a smaller lifetime, and 
hence a smaller efficiency. Second, a larger chargino mass leads to 
smaller stau momenta, and to smaller decay lengths. 

rriQ = 1000 eV/c^: In this case, the map of efficiencies is mainly affected by the mo- 
mentum of the stau, because the method used to identify heavy 
stable particles relies on the lack of Cherenkov radiation in DEL- 
PHI'S RICH detectors. To remove SM backgrounds, low momen- 
tum particles are removed, thus reducing the efficiency for higher 
chargino masses, especially in the region of small Am. 

4.4 Interpretation 

4.4.1 Neutralino pair production 

Given the aforementioned limits for the production cross- sect ion, some sectors of the 
(m^o , m^-) space can be excluded. In order to achieve the maximum sensitivity, the results 
from two other analyses are taken into account. The first is the search for slepton pair 
production in the context of gravity mediated SUSY breaking models |T^. In the case 



where the MSUGRA is massless, the kinematics corresponds to the case of / decaying 
into a lepton and a gravitino. The second is the search for lightest neutralino pair 
production in the region of the mass space where Xi is the NLSP (the region above 
the diagonal line in fig. |TD|, i.e. rrif > m^o). Within this zone, the neutralino decays into 
a gravitino and a photon. 

As an illustration, fig. [1^ presents the 95% C.L. excluded areas for rriQ < 1 eV/c^ in 
the m^o vs. mj^ plane for the co-NLSP case. The positive-slope dashed area is excluded 
by this analysis. The resulting 95% C.L. lower limit for the mass of the lightest neutralino 
is 82.5 GeV/c^. The negative-slope dashed area is excluded by the analysis searching for 
neutralino pair production followed by the decay Xi ~^ G'y. The point-hatched area is 



excluded by the direct search for slepton pair production within MSUGRA scenarios |12 . 
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Table ^ shows the 95 % CL. lower limits on the mass of the lightest neutralino within 
the two scenarios for different n. 

4.4.2 Slepton pair production 

The fi pair production cross-section depends on the mixing in the stau sector. There- 
fore, in order to put limits to the fi mass the mixing angle has to be fixed. The results 
presented here corresponds to two cases within the fi-NLSP scenario. In the first case, it 
is assumed that there is no mixing between the fn and fi- Thus, fi is a pure right-handed 
state (figure The second case (figure |ll]-b), corresponds to a mixing angle which 

gives the minimum fi pair production cross-section, while at the same time holds m?^ > 
0. Since this angle is close to the region where the coupling to the almost vanishes, no 
limit can be inferred from the LEPl measurements, and the search was extended down 
to stau masses around 2 GeV/c^. Stable or long lived particles with masses down to 
2 GeV/c^ are excluded by the search for heavy stable and long-lived particles in DEL- 
PHI [|TT|. Stable or long lived particles with masses below 2 GeV/c^ are excluded by the 
JADE collaboration RTI]. It is assumed that a stau with lower mass than a tau is stable 
or very long lived. In the case of short lived staus {rriQ < 0.03 eV/c^), a narrow band 
at rrir < nif^ < 2 GeV/c^ is not excluded. Above 2 GeV/c^ the results from the impact 
parameter and secondary vertex analyses are used for exclusion purposes. 

The impact parameter and secondary vertex analyses allow for the exclusion of fi (f/j) 
with a mass below 80 GeV/c^ for gravitino masses between 10 and 310 eV/c^ (8 and 
380 eV/c^) at 95% C.L.. For below a few eV/c^, < 73 GeV/c^ were excluded by 
the search for fi in gravity mediated models [jl2|. Results from both searches were not 
combined because the impact parameter searches cover in excess the overlaping region. 
For tjIq larger than 1000 eV/c^ the limit is 87 GeV/c^, obtained after combining the 



results presented in this paper with those of the stable heavy lepton search [|ri[ . 

Within the sleptons co-NLSP scenario, the cross-section limits of figures ^b and -c 
were used to derive limits for /t^ (fig 0-a) and (fig. |12|-b) masses at 95% C.L.. The 
fiji- pair production cross-section is model independent since it only takes place through 
the exchange of a Z*^ or a 7 in the s-channel. The en- pair production cross-section, 
however, is a function of the GMSB parameters due to the exchange of a Xi the t- 
channel. Therefore, in order to put limits to the en mass, the aforementioned region of the 
GMSB parameter space was scanned, and for each selectron mass the smallest theoretical 
production cross-section was chosen for comparison with the experimental limits. 

The t-channel interference causes a bigger fraction of selectrons to be produced in 
the forward region. This results in a loss of efficiency in the vertex and stable slepton 
analyses, that was taken into acount for the calculation of the limits that are shown in 
fig. 0-b. 

Therefore, within the co-NLSP scenario, the impact parameter search and the sec- 
ondary vertex search allow for the exclusion of ftR masses below 80 GeV/c^ for gravitino 
masses between 8 and 450 eV/c^. In the case of cr, masses below 67 GeV/c^ for gravitino 
masses between 10 and 80 eV/c^ are excluded. 

Assuming mass degeneracy between the staus and smuons, (fig. ^2]-c), these searches 
exclude at 95% C.L. Ir masses below 84 GeV/c^ for G masses between 9 and 570 eV/c^. 
For very short lifetimes only JIr was considered since it is the best limit that can be 
achieved in absence of slepton combination. For G larger than 1000 eV/c^ the limit was 



87 GeV/c^, obtained from the stable heavy lepton search |TT|]. Ir masses below 35 GeV/c 
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are excluded from LEP 1 data |32|. In the case of //j degeneracy, this hmit improves to 
41 GeV/c^. 

4.4.3 Chargino pair production 

The hmits on chargino pair production cross-section were used to exclude areas within 
the (m-+,mfj plane in different domains of the gravitino mass @]. 



Figure |T3| shows the regions excluded at 95% CL in the (m-+,m^J plane. The positive- 
slope area is excluded for all gravitino masses. The negative-slope area is only excluded 
for rriQ > 100 eV/c^. The area below m^^ = 75.8 GeV/c^ is excluded by the direct search 
for stau pair production [12|. The area of Am < 0.3 GeV/c^ is not excluded because in 
this regions the charginos do not decay mainly in fi and Ut, but in W and G. Thus, if 
Am > 0.3 GeV/c^, limits at 91.8, 93.0 and 93.0 GeV/c^ can be set for = 1, 100 and 
1000 eV/c^ respectively. The limit at fnQ = 1 eV/c^ is also valid for smaller masses of the 
gravitino, because they lead to the same final state topologies. The same argument is 
true for m^, >1 keV/c^. The chargino mass limit decreases with decreasing m^^ because 
in scenarios with gravitino LSP, small stau masses correspond to small sneutrino masses 
(both are proportional to A), and hence to smaller production cross-sections due to the 
destructive interference between the s- and t-channels. It should be noticed that within 
the parameter space that concerns this work, the lightest chargino is at least 40% heavier 
than the lightest neutralino. Thus, for small gravitino masses the search for neutralinos 
implies a lower limit on the lightest chargino of 120 GeV/c^. Neutralinos are not directly 
searched for in heavier gravitino mass regions and therefore the limit of 93 GeV/c^ remains 
valid. 



4.4.4 Limits on the GMSB parameter space 

Finally, all these results can be combined to produce exclusion plots within the 
(tan/5,A) space. As an example, fig. shows the zones excluded for =1 to 4 for 



iTT-G — 1 eV/c , which corresponds to the NLSP decaying at the main vertex. The shaded 
areas are excluded. The areas below the dashed lines contain points of the GMSB param- 
eter space with x^-NLSP. The areas to the right (above for n = 1) of the dashed-dotted 
lines contain points of the GMSB parameter space where sleptons are the NLSP. It can 
be seen that the region of slepton-NLSP increases with n. The contrary occurs to the 
region of neutralino-NLSP. A limit can be set for the variable A at 16.5 TeV. 



5 Summary 

Lightest neutralino-, slepton- and chargino pair production were searched for in the 
context of light gravitino scenarios. Two scenarios were explored: the fi NLSP and 
the Ifi co-NLSP scenarios. No evidence for signal production was found. Hence, the 
DELPHI collaboration sets lower limits at 95% C.L. for the mass of the Xi 82 GeV/c^ 
if niQ < 1 eV/c^, for the mass of the fi at 73 GeV/c^, the Ir at 79 GeV/c^, and the 
lightest chargino at 93 GeV/c^ for all rriQ. A limit is also set on the variable A at 16.5 
TeV if < 1 eV/c^. 
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Cut 


77 


//7 


4-fermion 


Total MC 


Data 


Signal 


pre-selection 


1134 ±37 


413 ± 14 


385 ± 18 


1933 ± 43 


1791 


59.5% 


1 


28 ±3 


330 ± 11 


363 ± 18 


721.3 ±21 


706 


51.2% 


2 


4.3 ± 1.2 


11.1 ± 1.5 


12.5 ±2.8 


28.0 ±3.4 


24 


42.1% 


3 





0.07 =b 0.07 


1.13 ±0.3 


1.2 ±0.3 


2 


32.8% 



Table 1: Number of events remaining in the data and simulated samples at ^/s — 189 GeV 
after the various stages of the selection procedure described in the search for neutralinos 
decaying into slepton and lepton. The signal efficiencies corresponds to m^o — 87 GeV/c^ 
and mi = 75 GeV/c^. 



Cut 


77 


ffl 


4-fermion 


Total MC 


Data 


Signal 


pre-selection 


1134 ±37 


413 ± 14 


385 ± 18 


1933 ± 43 


1791 


57.7% 


1 


66.5 ±7 


376.4 ± 13.0 


331.0 ± 12.1 


468.5 ± 19.1 


404 


54.5% 


2 


6.7 ± 1.8 


9.5 ± 1.3 


10.6 ±0.9 


26.8 ±2.4 


23 


44.9% 


3 





0.07 ±0.07 


1.09 ±0.3 


1.16 ±0.3 


1 


37.3% 



Table 2: Number of events remaining in the data and simulated samples at ^/s = 189 GeV 
after the various stages of the selection procedure described in the search for neutralinos 
decaying into stau and tau. The signal efficiencies corresponds to m^o — 87 GeV/c^ and 
rrif, = 75 GeV/c^. 



Observed events 


1 


Total background 


1 42+^-''^ 


Z*/^^{ll){rn) 


^•^'J-O.Ol 


4-fermion (except 77) 


0.45±0.16 


77 — > t'^t" 


'-'■'^-0.32 



Table 3: Number of observed events at y/s = 189 GeV, together with the total number of 
expected background events and the expected numbers from the individual background 
sources, for both large impact parameter and secondary vertex searches combined. 
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Observed events 


4 


Total background 




Z77^(rr)(n7) 


1 qq+U.46 
^^•'-''-'-0.35 


77 — i> r+r^ 


'-'•"-'--0.38 


WW 


z.jz_Q 23 


zz 





Table 4: Expected simulated SM background events and selected data events at 189 GeV 
centre-of-mass energy for the small impact parameter search. 



Sample 


Efficiencies (%) 


Main backgrounds 


Expected b.g. 


Observed events 


uiq = 1 eV/c^ 


24 - 36 


WW , 77 


38.9±4.9 


36 


= 100 eV/c^ 


28 - 50 


77 


2.1±0.9 


1 


rric = 1000 eV/c^ 


- 63 




1.7±0.3 


1 



Table 5: Range of efficiencies for the different sets of chargino signals described in sec- 
tion 1^, main sources of background, expected background and observed data events for 
the different analyses. 



n 


co-NLSP( GeV/c' ) 


n-NLSP ( GeV/c' ) 


1 


83.0 


82.5 


2 


85.0 


85.0 


3 


86.0 


86.0 


4 


87.0 


87.0 



Table 6: The 95% CL. lower limits on m^o within the Ti-NLSP and co-NLSP scenarios 
for different n. 
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Figure 1: Two examples of kinematic differences between the fi-NLSP and co-NLSP 
scenarios. Fig. (a) shows the distribution of the missing mass normahzed to the cen- 
tre-of-mass energy {rrimiss/y/s) for simulated sets with same mass neutralinos and same 
mass fi and slepton. Fig. (b) shows the number of charged tracks per event for the 
same two sets of simulated signals. Histograms with positive-slope shading show a set of 
m^o = 82 GeV/c^ and ruf^ = 80 GeV/c^. Histograms with negative-slope shading show 
a set of m^o = 82 GeV/c^ and mj^ = 80 GeV/c^. 
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Figure 2: (a) Normalized missing mass and (b) momentum of the leading charged particle, 
for data (dots), Standard Model simulation (cross-hatched histogram) and one of the 
simulated signals with cross-section with arbitrary normalization (blank histogram) after 
preselection. The arrows indicate selection criteria imposed as explained in the text. 
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Acoplanarity (") 

Figure 3: Acoplanarity of data (dots), Standard Model background simulation 
(cross-hatched histogram) and one of the simulated signals with cross-section not to 
scale (blank histogram), after the cut to remove 77 events. The arrow indicates selection 
criterion imposed as explained in the text. 
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Figure 4: (a) Minimum angle of a jet to the beam, (b) maximum of angular broadness of 
the two jets at the 2-jet level and (c) minimum angle between jets at the 4-jet level, for 
data (dots). Standard Model background simulation (cross-hatched histogram) and one 
of the simulated signals with cross-section not to scale (blank histogram), after the cut 
to remove //(7) and 4-fermion events. The arrows indicate selection criteria imposed as 
explained in the text. 
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Angle with shower (") 



Figure 5: (a) Angle between the directions defined by the hadronic vertex and the re- 
constructed vertex, (b) angle between the tracks of the kink, and (c) angle between the 
electromagnetic shower and the direction defined by the difference between the momenta 
of the fi and its associated r, defined at the crossing point for real data (dots), ex- 
pected Standard Model background (cross-hatched histogram) and simulated signal for 
rrif^ = 60 GeV/c^ decaying with a mean distance of 50 cm (blank histogram). Events 
that do not have hadronic interactions are not included in (a), and events without elec- 
tromagnetic showers are not included in (b). The arrows indicate the selection criteria 
imposed. 
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Figure 6: Missing transverse momentum (a) and visible energy (b), for real data 
(dots), expected Standard Model background (shaded histogram) and simulated signal 
for = 85 GeV/c^ and m-+ = 69 GeV/c^ decaying with a mean distance of 50 cm 
(blank histogram). The arrows indicate selection criteria imposed as explained in 
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Figure 7: 95% C.L. upper limit of the x? pair production cross-section (in femtobarn) at 
^/s = 189 GcV (a) after combining the results of the searches from ^/s = 161 up to 189 
GeV, as a function of m^o and m^^ for the case n = 3 in the fi-NLSP scenario, where 
n is the number of messenger generations and (b) using data at y/s = 189 GeV, as a 
function of m^o and nii^ in the co-NLSP scenario. The diagonal and vertical lines show 
respectively the limits m^o = mf + rrir and m^o — y/s/2. 
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Figure 8: 95% C.L. upper limit of the e'^e~ — > fifi (a), e'^e~ — > /iR/ii? (b), e'^e~ — > e^eR 
(c), and e~^e~ IrIr (d) production cross-sections, at y^=189 GeV after combining 
the resuhs of the searches at = 130—189 GeV. Results are shown in the (m^,m;-) 
plane. Searches for events containing charged particle tracks with small impact param- 
eter, large impact parameter, secondary vertices and the search for heavy stable leptons 
are combined. 
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Figure 9: Limits in picobarn on the lightest chargino pair production cross-section at 
95% CL. Limits are shown as functions of m-+ and ruf^ for (a) uIq — 1 eV/c^, (b) 
ma - 100 eV/c2 and (c) = 1000 eV/c^. 
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Figure 10: Areas excluded at 95% C.L. with m^^ < 1 eV/c^ in the m^o vs. m^-^ plane for 
n =1 to 4. The positive-slope dashed area is excluded by this analysis. The negative-slope 
dashed area is excluded by the search for Xi ~^ iG, and the point-hatched area by the 
direct search for slepton pair production in the MSUGRA framework. 



25 



p I I I I I I I I I I I I I I I 1 1 I I I I I I I I I I I I I I I I I 



:^ 90 

<u 



80 
70 
60 
50 
40 



MSUGRA ^// 




Stable 
Heavy 
Lepton n 



10 



10^ 10^ 

2 

mg (eV/c ) 



90 

^ 80 
O 

^^70 
S 60 
50 
40 
30 
20 
10 




hi 1 1 1 1 1 I I I 1 1 1 1 1 1 I I I 1 1 1 1 1 1 I I I 1 1 1 1 1 1 I I I 1 1 1 1 1 1 I I I 1 1 1 1 1 1 I I I 1 1 1 1 1 



1 



(b) 



MSUGRA 




10 



10 



-2 



10 



10 



stable 
Heavy 
Leptoii 



10 ' 10 ' 

mg (eV/c ) 



Figure 11: Exclusion region in the {mQ,mfj^) plane (a) and {mQ^nif^) plane (b) at 
95% C.L. for the present analysis combined with the stable heavy lepton search and 
the search for fi within MSUGRA models, using all LEP-2 data up to 189 GeV. The 
positive-slope hatched area shows the region excluded by the impact parameter and 
secondary vertex searches. The dashed line shows the expected limits. The area below 
the horizontal line in (b) is excluded by the JADE collaboration. A narrow band at 
ttIt < nif < 2 GeV/c^ and uiq < 0.03 eV/c^ is not excluded. 
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Figure 12: Exclusion regions in the (m(5,m^^) (a), (m^jjme^) (b) and {mQ,mj^) (c) 
planes at 95% CL. for the present analyses combined with the stable heavy lepton search 
and the search for in gravity mediated models, using all LEP2 data up to 189 GeV. 
The positive-slope hatched area shows the region excluded by the combination of the 
impact parameter and secondary vertex searches. The dashed line shows the expected 
limits. 
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Figure 13: Areas excluded at 95% CL in the (m^+,mfj plane. The positive-slope area is 
excluded for all m^. The negative-slope area is excluded only for > 100 eV/c^. The 
grey area is excluded by the direct search for stau pair production . 
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Figure 14: Shaded areas in the (tan (3, A) plane are excluded at 95% CL. The areas below 
the dashed lines contain points of the GMSB parameter space with x^'NLSP. The areas 
to the right (above for n = 1) of the dashed-dotted lines contain points of the GMSB 
parameter space were sleptons are the NLSP. 



